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C^ , Abstract: This paper introduces a new distribution to improve tail risk modelling. Based 

f^ , on the classical normal distribution, we define a new distribution by a series of heat equations. 

^^ I Then, we use market data to verify our model. 
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^^ ; 1 Introduction 

r^ I According to Espen Gaarder Haug(2007), Wesley C. Mitchell conducted the first empirical 
study of fat-tailed (high-peaked) distributions in price data as early as 1915. Following Man- 






delbrot's famous 1962/63 paper on fat-tails, there have been an abundant research focusing 
on the non-explained empirical fact. Time-varying volatility, standard returns, mixture model, 
jump-diffusion model, stochastic volatility, implied distributions are all important tools for risk 
measurement and management. However, the 2007-08 financial crisis exposed the deficiencies 
^ ' in existing risk models and reinforced the importance of methodology improvements. More 

Q> ■ detialsee [T], [2], [2]. 

In this paper, we consider a series of heat equations which relate with normal distribution. 
■ ■ Assuming the volatility of random varible dependent on the value of random varible, we define 

f^ ■ a new distribution-random limit normal distribution. 
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The remainder of the paper followed as: Section 2 introduces a new distribution-random 
limit normal distribution, and gives two usefull parameter models. Section 3 use the random 
limit normal distribution to fit history data. Some technique proof is given in the Section 4. 

2 A New Distribution 

2.1 Random Limit Normal Distribution 

Alternatively, we propose a definition of normal distribution: 



'Institute for Economics and Institute of Financial Studies, Shandong University, Jinan, Shandong 250100, 
PR China(agcaelyn@gmail.com). 

'School of mathematics, Shandong University, Jinan, Shandong 250100, PR China(yangsz@mail. sdu.edu. en). 



Definition 2.1 A -dimensional random vector X = [Xi,--- ,X(i) on a sublinear expectation 
space {Q,T,P) is called (centralized) normal distribution, if 



aX + bX = Va^ + b'^X, Va, 6 > 0, 

where X is an independent copy of X, whicli means X is independent of X , and X ^ X . 
When d = 1, wc have X = N{0, cr^), where a^ = E[X^]. 

Let u{t,x) — E[(f{x + \/tX)], then u{t,x) satisfies the following heat equation: 

dtu{t, x) - -a^dl^u{t, x) = 0, m(0, x) = ip{x), x e R. (2.1) 

Set (p{x) = l{j.<j,-j., we have 

E[^{X)] =P(X <y)= uy{l,0) = -=L= r exp(-f^)dx. (2.2) 

So the normal distribution of X is 

Fx(y)-u,(l,0), yei?. 

Following the above definition of normal distribution, we give the definition of the random 
limit normal function. 

Definition 2.2 For a given X on a probability space {fl,T,P), with function h{-) : R — > 
(0,+cx)), and the heat equation 

dtu{t,x)- -a^h^(y)dl^u{t,x) =0, u{0,y) = l^^^yy, x e R. 

We define the random limit normal function Fx{-) of X as 

Fx{y) - uy{l, 0) = , }^ ^^^ f exp(-— ^^^)dx, (2.3) 

v/27r(cr/i(y))2 7_oo 2(cr/i(y))2 

where EX ^Q, EX^ = a^ . 

Remark 2.3 In general, the random limit normal function is not a distribution. However, in 
the following, we will show two cases when it become a distribution . 

Theorem 2.4 For a given X on a probability space {n,T,P), EX = 0, EX^ = a^ . For 
function h{-) : R — 5- (0, +oo), // any of the following case is right, Fx{-) is a distribution. 

case 1: h{-) is continuous increasing on (— oo,0), and continuous decreasing on [0, +oo); 

case 2: Set h{y) = Ky + c, c> 0, if y < 0, K < 0, if y > 0, K > 0. 

Proof, case 1: If y < 0, for a given small (5 > 0, set y + 5 < 0. Then, by the definition of Fx{-), 
we have 



1 
y/27r(ah{y+5)y 



p^(y + s)= ./2.(.l+.-))^ /-^'^"p(- 2(..g+^))- )^"- 



Note that, (Th{y + S) > ah{y). By the comparing of two normal distribution, we have 

Fx{v + 5)>Fx{y). 
Similar, if y > 0, we aslo have 

Fx{y + 5)>Fx{v). 

So Fx{-) is a increasing function on R, and when y — > +00, Fx{y) — > 1- 
The proof of case 2 is more technique, we will show it in the Appendix. ■ 

Remark 2.5 Note that, we don't need the condition EX = m the proof of Theorem \2.4\ 
From now, we use the function h in the case 2, i.e. h{y) = K ■ y + c. For reader convenience, 
we denote N(u,a'^, K,c) as the random limit normal distribution of X. 

3 Two Example 

3.1 Test of Hypothesis 

Following the classical work of K. Pearson, we conduct goodness of fit test for random limit 
normal distribution. 

When use the method to group the data, the number of groups is m, while n is the sample 
size. We suppose there are mutually disjoint intervals /i, • • • ,Im- 

m Yi 

A=E-fe-P»)' (3.1) 

i=iPi 

where qi is the frequency of the i tli group, 1 < i < m, andpi is the corresponding probability 
of the fitting distribution. 

K. Pearson proves the following theorem: 

Theorem 3.1 // the theoretical distribution is right, when the size of sample n — > 00, the 
limit distribution of statistic A is the y^ distribution in {Q,J^,P) with k — 1 degrees of freedom. 

If the parameters in theoretical distribution are unknown, according to the theorem of 
Fisher, we take the limit distribution of statistic A as the x^ distribution in (i7,J^, P) with 
fc — 1 — r degrees of freedom. And r is the number of unknown parameters. 

3.2 Afitting Example 1 

For distributions N{iJL,a'^) and N{iJL,a^ ,K,c), we estimate /i and a by classical method. The 
fitting result of normal distribution and random limit normal distribution is given. We use the 
test index 5 and 5'^"^ to compare normal distribution and random limit normal distribution, 5 
and (5^^" follow as: 

m r, 

5=Y.-{<i^-P^? (3.2) 

i=lPi 



i=lPi 



(3.3) 



Where qi is the frequency of the i th group, 1 < i < tti, m is the number of group, n is the 
number of data, pi is the corresponding probability of the fitting normal distribution, p^*^" is 
the corresponding value of the fitting random limit normal distribution. 

Using N{fi,a^) to fit the data of S@P 500 2000-2012 Log Return, and the parameters is 

The style of data S@P 500 

Length of data 3311 

fi 1.70610051464811e-005 

a 0.0134453931516791 

S 1129.40451908248 

and the fitting figure is 




Using N{^,a'^,K,c) to fit the data of S@P 500 2000-2012 Log Return, and the parameters 



IS 



The style of data 


S@P 500 


Length of data 


3311 


^^ 


1.70610051464811e-005 


a 


0.0134453931516791 


K 


-24(y<0),24(y>0) 


c 


0.5 


jccn 


349.545332843509 



and the fitting figure is 




Using N{^i,a'^) to fit the data of S@P 500 1950-2012 Log Return, the parameters is 

The style of data S@P 500 

Length of data 15852 

fi 1.70610051464811e-005 

a 0.0134453931516791 

6 4830.10362738027 
and the fitting figure is 




-0.04 -0.02 0.02 0.04 0.06 0.08 0.1 



Using iV(^, (t'^,K, c) to fit the data of S@P 500 1950-2012 Log Return, tlie parameters is 



The style of data 
Length of data 

/^ 
a 
K 



S@P 500 

15852 

1.70610051464811e~005 

0.0134453931516791 

-30(2; <0), 30(2; >0) 

0.56 

1201.24543468315 



and the fitting figure is 
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3.3 Afitting Example 2 

Using 7V(/i, 0-2) to fit the data of MSFT 1986-2012 Log Return, the parameters is 

The style of data MSFT 

Length of data 15852 

/I 1.70610051464811e-005 

a 0.0134453931516791 

5 2388.80170177018 



and the fitting figure is 




Using N{^, CT^, K, c) to fit the data of MSFT 1986-2012 Log Return, the parameters is 



The style of data 
Length of data 

a 
K 



MSFT 

15852 

1.70610051464811e-005 

0.0134453931516791 

-28.3(2; <0),28.3(j/>0) 

0.24 

1201.24543468315 



and the fitting figure is 




Recent research is strongly motivated by the financial meltdown to develop models and 
methods to better predict extreme events. In this paper, we propose a new distribution: con- 
tinuous change normal distribution which includes the normal distribution as a special case 



and takes into account of main statistically relevant stylized facts of asset returns data: fat-tail, 
high-peak. 

4 Appendix 

The proof of case 2 in Theorem 12. 4t 

Proof, case 2: For convenience, we set a — 1. 

when y e (— oo, 0], K < 0. In order to prove Fx{-) is a distribution, we need to show that 
Fx{-) is a continuous and increasing to I's function. By the definition of Fx{-), we have 

Considering the right derivative function of Fx{-), and denote as fx{')- 



(4.1) 



i/27r(Kjj+c)2 "^^P^ ■2(Ky+cy^- 



Then, we have 



/x(0) 



K I K 



Set 



^27r(A'y+c)2 2(Ky+c) 2(Kv+c) 

_ 1 

^ ^27r(Xy+c)2 

> 0. 



Set ?7i(y) = K^+c ^ then m G (-^,0]. By parameter transformation, we have 

~gx{fn) = exp(-il^) +/:"^if[exp(-4) • (z^ _ i)]dx. 

Note that, ttt, is a increasing to O's function of y. Also we have 5x(0) > 0. The derivative 
function oigx{-) is 

-g'xim) = - exp(-^) • m + i^[exp(-^) • {m^ - 1)] 
= exp(-2|^) • (Km2 - m ~ X). 

For m £ {^,0], so g'xi') > 0. Then 5x(') is increasing function of m, and also gx{') and 
/x(') is a increasing function of y. 

By the equation (|4.ip . fx{y) — >■ 0, when y — > —00. So fx{y) > 0, y € (— oo,0], and Fx{-) 
is a increasing function on (—00, 0]. 



Similar, we could prove Fx{-) is a increasing function on R, and when y — ;■ +00, Fx{y) — > 
1. 

This complete the proof. ■ 
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